AD-A237  462 


fATION  PAGE 


form  Approved 
OM8  No  0704  9>*fl 


1b  RESTRICTIVE  MARKINGS 


?b  DECLASSIFICATION  I  DOWNGRADING  SCHED 


a  PERfORMING  ORGANIZATION  REPORT  NUMB£R(S) 


3  DISTRIBUTION/ AVAILABILITY  Of  REPORT 

Approved  for  public  release  and  sale;  its 
distribution  is  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6b  OffICE  SYMBOL 
(If  applicable) 


8b  OffICE  SYMBOL 
(If  applicable) 


Technical  Report  No.  101 


6a  NAME  Of  PERfORMING  ORGANIZATION 
Purdue  University 
Department  of  Chemistry 


6c.  ADDRESS  (City.  State,  and  IIP  Code) 
Purdue  University 
Department  of  Chemistry 
West  Lafayette.  IN  47907 


Da  NAME  Of  FUNDING /SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 


8c  ADDRf  SS  (City,  State,  and  ZIP  Code) 

800  N.  Quincy  Street 
Arlington,  VA  22217 


il  HUE  (include  Security  Claudication) 

Solvent  Dynamical  Effects  in  Electron  Transfer:  The  Solvent  Inertial  Limit  and  the 
Predicted  Influence  of  Quantum  Effects 


1 2  personal  AUTHOR(S) 

G.E.  McManis,  A.  Gochev,  and  M.J.  Weaver 


13a  TYPE  OF  REPORT  113b  TIME  COVERED 

Technical  |  from  to 


7a  NAME  OF  MONITORING  ORGANIZATION 
Division  of  Sponsored  Programs 

Purdue  Research  Foundation 


7b.  ADDRESS  (Ciry.  State,  and  ZIP  Code) 
Purdue  University 
West  Lafayette,  IN  47907 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMPf  R 

Contract  No.  N00014-91-J-1409 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  1  PROJECT 

ELEMENT  NO  I  NO 


14.  DATE  OF  REPORT  (Year,  Month.  Day)  15  PAGE  CuUNl 

May  31,  1991 


I  7. _ COSATI  CODES _ I  18  SUBJECT  TERMS  (Continue  on  reverse  >f  necessary  and  identify  by  block  number) 

ri£LD  1  group  |  subgroup  I  analytical  continuum  formulations  for  the  solvent  inertial 

frequency,  nuclear  tunneling,  rate-solvent  friction 


19  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Analytical  continuum  formulations  for  the  solvent  inertial  frequency,  constituting  cr.e 
anticipated  zero-friction  limit  for  adiabatic  barrier  crossing  in  solvent-controlled 
electron-transfer  processes,  are  derived  and  discussed.  The  role  of  solvent  inertia  in 
solvent  dynamical  effects  is  discussed  with  emphasis  on  the  likely  modifications  brought 
about  by  nuclear  tunneling.  Approximate  formulations  suitable  for  assessing  the  latter 
correction  in  the  presence  as  well  as  absence  of  solvent  friction  are  outlined  and  compared. 
Numerical  calculations  are  provided  that  illustrate  the  partly  compensatory  influence  on  the 
rate-solvent  friction  dependence  resulting  from  nuclear  tunneling  together  with  reaction 
nonadiabaticity .  Such  combined  quantum  effects  are  anticipated  typically  to  mask  the 
clearcut  emergence  of  solvent  inertial  effects  on  electron-transfer  reaction  dynamics  in 
common  low-friction  media. 


Z0  DISIRIBUTION/AVAILABILTTY  OF  ABSTRACT 
□  UNCLASSIFIED/UNUMITED  □  SAME  AS  RPT  Q  OTIC  USERS 


}2*  NAME  OF  RESPONSIBLE  INDIVIDUAL 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


22b  TELEPHONE  (include  Area  Code)  22c  OMiCE  SyMR<h 


DO  form  1473,  JUN  86 


Previous  editiont  are  obsolete 


SECURITY  CLASSiF If  ATIQN  Ot  this  rsi  > 


OFFICE  OF  NAVAL  RESEARCH 


Contract  No.  N00014-91-J-1409 
Technical  Report  No.  101 


Solvent  Dynamical  Effects  in  Electron  Transfer: 

The  Solvent  Inertial  Limit  and  the  Predicted  Influence  of  Quantum  Effects 


by 


G.E.  McManis,  A.  Gochev,  and  M.J.  Weaver 
Prepared  for  Publication 
in 


Chemical  Physics 


Purdue  University 
Department  of  Chemistry 
West  Lafayette,  Indiana  47907 


May  1991 


■  fc-'S  QKuiI 

;  vsic 

j  Uaanoo'vitied 
•  Justification 


Distribution/ 


- . -  r  _ .  . 

Availability  Codee 

Diet 

M 

Avail 

Spa« 

• 

13 

Reproduction  in  whole,  or  in  part,  is  permitted  for  any  purpose  of  the  United 
States  Government. 

*  This  document  has  been  approved  for  public  release  and  sale:  its  distribution 
is  unlimited. 


Solvent  Dynamical  Effects  in  Electron  Transfer: 
The  Solvent  Inertial  Limit  and  the  Predicted 
Influence  of  Quantum  Effects^ 


George  E.  McManis,  Alexander  Gochev, 
and  Michael  J.  Weaver 

Department  of  Chemistry,  Purdue  University 
West  Lafayette,  Indiana  47907  (USA) 


Chemical  Physics 

(Special  Issue:  "Friction  in  Liquid-Phase  Reactions") 

Submitted  May  30,  1990 
Revised  October  15,  1990 

^Abstracted  in  part  from  Ph.D.  thesip  of  G.  E.  McManis,  Purdue  University, 
1989 


abstract 


Analytical  continuum  formulations  for  the  solvent  inertial  frequency, 
constituting  the  anticipated  zero- friction  limit  for  adiabatic  barrier 
crossing  in  solvent-controlled  electron- transfer  processes,  are  derived  and 
discussed.  The  role  of  solvent  inertia  in  solvent  dynamical  effects  is 
discussed  with  emphasis  on  the  likely  modifications  brought  about  by  nuclear 
tunneling.  Approximate  formulations  suitable  for  assessing  the  latter 
correction  in  the  presence  as  well  as  absence  of  solvent  friction  are  outlined 
and  compared.  Numerical  calculations  are  provided  that  illustrate  the  partly 
compensatory  influence  on  the  rate -solvent  friction  dependence  resulting 
from  nuclear  tunneling  together  with  reaction  nonadiabaticity .  Such  combined 
quantum  effects  are  anticipated  typically  to  mask  the  clearcut  emergence  of 
solvent  inertial  effects  on  electron- transfer  reaction  dynamics  in  common 
low-friction  media. 
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I.  INTRODUCTION 

The  realization  that  the  dynamical  properties  of  the  solvent  can  play 

an  important  role  in  the  kinetics  of  electron  transfer,  as  in  other  condensed- 

phase  chemical  processes,  has  spawned  an  impressive  range  of  recent  studies 

of  both  an  experimental  and  theoretical  nature  (see  refs.  1  and  2  for  recent 

overviews) .  Inquiries  of  the  former  type  have  encompassed  examinations  of 

activated  electron- transfer  (ET)  reactions  (i.e.  featuring  significant  free- 

energy  barriers),^"  as  well  as  real-time  dynamics  of  polar  solvation  as  studied 

by  ultrafast  laser  techniques,  especially  time -dependent  fluorescence  Stokes 

2  3 

shift  (TDFS)  measurements.  ’  These  considerations  have  led  to  the  conclusion 
that  in  many  solvents,  at  least  for  near-ambient  conditions,  the  nuclear 
reaction  dynamics  are  determined  primarily  by  overdamped  solvent  relaxation 
as  characterized  in  the  simplest  case  by  the  longitudinal  relaxation  time, 
rL .  In  other  words,  the  common  occurrence  of  solvent  friction  depresses 
the  net  frequency  of  adiabatic  barrier  crossing  below  the  value,  w,/2*  (s'1), 
reflecting  the  inertia  of  individual  solvent  dipoles.  The  latter  frequency 
is  that  expected  on  the  basis  of  transition- state  theory  (TST)  treatments. 

Nevertheless,  it  is  anticipated  that  in  some  circumstances  the  solvent 
friction  will  be  sufficiently  small  so  that  tl'x>  ,  whereupon  the  rate  of 
barrier  crossing  can  be  determined  partly  or  wholly  by  the  longitudinal 
inertial  polarization  dynamics  as  reflected  in  w0 ,  rather  than  by  the  over¬ 
damped  frequency  r  ^  1  .  In  principle,  inertial-limiting  effects  could  be 
discerned  by  non- exponential  TDFS  decays  at  short  times^  and  also  from  rL  - 
independent  barrier-crossing  frequencies  for  suitably  adiabatic  ET  processes.5 
In  contrast  to  overdamped  motion,  however,  no  clearcut  experimental 
observation  of  such  inertial  effect?  ir.  polar  solvents  has  apparent 1 y  been 
reported.  Such  observations  are  hampered  by  the  high  frequencies  anticipated 
for  solvent  dipole  inertia,  typically  u>0  >  5  x  1012  s'1  for  typical  dipolar 
solvents.  Since  thereby  u>0  >  k^T/ft  (where  kg  is  the  Boltzmann  constant)  at 
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ambient  temperatures,  quantum  effects  may  commonly  be  anticipated,  specifi¬ 
cally  involving  nuclear  tunneling  so  to  enhance  the  barrier-crossing  rate 
above  the  expectations  of  classical  dynamics.  An  additional  major  limitation 
to  the  emergence  of  solvent  inertial  effects  is  that  the  maintenance  of  reac¬ 
tion  adiabaticity ,  whereupon  the  full  effect  of  nuclear  dynamics  upon  the 

ET  rates  will  be  felt  even  at  such  rapid  frequencies,  requires  large  degrees 

6  7 

of  donor -acceptor  electronic  coupling.  ’  Conversely,  diminishing  values 
of  the  electronic  transmission  coefficient,  <c#1,  are  often  expected  as  the 
friction  is  decreased  by  altering  the  solvent,**  so  that  such  reaction  non- 
adiabaticity  can  often  curtail  the  influence  of  solvent  inertia  on  the  bar¬ 
rier-crossing  dynamics  under  these  conditions.  Consequently,  then,  the  occur¬ 
rence  of  inertial  limiting  effects  on  solvent -dependent  ET  reaction  rates 
may  be  obscured  by  the  increasingly  accelerating  and  retarding  influences 
of  nuclear  tunneling  and  reaction  nonadiabaticity ,  respectively,  as  the  nu¬ 
clear  dynamics  become  more  facile. 

The  central  issue  addressed  herein  is  the  numerical  assessment  of  the 
likely  role  and  importance  of  solvent  inertia  on  ET  reaction  rates  given 
the  inevitable  additional  presence  of  such  quantum  effects.  The  organization 
of  the  paper  is  as  follows.  We  first  discuss  specific  analytic  formulations 
for  the  inertial  limiting  frequency  and  related  quantities,  so  to  characterize 
the  magnitudes  and  physical  origins  of  such  effects,  and  to  provide  approxi¬ 
mate  estimates  of  «o  .  The  latter  is  of  practical  interest  given  the  absence 
of  direct  experimental  information  on  inertial  frequencies.  We  then  outline 
theoretical  formulations  that  enable  the  nuclear  quantum  corrections  to  be 
made  in  the  presence  as  well  as  absence  of  solvent  friction.  This  is  followed 
by  representative  numerical  calculations  so  to  illustrate  the  likely  experi¬ 
mental  consequences  of  coupled  inertial-quantal  effects  upon  ET  exchange 
rates  in  low- friction  media.  Of  particular  interest  here  is  the  manner  in 
which  these  factors  are  liable  to  influence  rate-solvent  friction  dependencies 
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since  this  issue  is  of  specific  relevance  to  recent  (and  ongoing)  experimental 
studies.^’ ^  While  the  emphasis  here  is  on  the  practical  consequences  of  such 
effects,  some  explanative  discussion  of  the  underlying  physical  origins  is 
provided  since  we  found  the  theoretical  literature  on  this  topic  often  to 
be  as  esoteric  as  it  is  extensive. 


II.  SOLVENT  INERTIAL  EFFECTS  ON  THE  ELECTRON -TRANSFER 


CONSTANT 


1.  Inertial  Limiting  Frequency 


The  unimolecular  rate  constant  for  an  electron- transfer  reaction  (such 

as  within  a  given  internuclear  reactant  geometry  for  a  bimolecular  outer- 
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sphere  reaction  )  is  conveniently  expressed  as 

k.t  “  rn  exP  (AG'AbT)  (la) 

-  rn  *elkA  <lb> 

where  t/n  is  the  nuclear  frequency  factor  and  AG*  is  the  activation  free  ener¬ 
gy.  The  terms  Fn  and  »c#1  provide  quantum  corrections  to  the  "classical" 
adiabatic  rate  constant  kA  arising  from  nuclear  tunneling  and  electronic 
nonadiabaticity ,  respectively.  For  a  given  reaction,  most  simply  for  a  fixed 
AG*  and  for  a  given  donor -acceptor  electronic  coupling  (as  gauged  by  the 
electronic  matrix  coupling  element  H1Z),  Tn  will  tend  to  increase  and  «el 
to  decrease  for  progressively  larger  values  as  engendered  most  readily 
by  altering  the  solvent.  For  suitably  small  ,  on  the  other  hand,  both  Tn 
and  »c#1  approach  unity  so  that  k>t  =>  kA  . 

In  high- friction  ("overdamped")  solvents,  the  adiabatic  barrier-crossing 
frequency  i/n  is  characterized  in  the  continuum  limit  by  r"1,  this  frequency 
can  fall  substantially  below  the  TST-limit  value,  «o/2x.  Nevertheless,  theo¬ 
retical  treatments  suggest  that  vn  is  not  greatly  different  from  wo/2 n  in 
some  media,  so  that  significant  solvent  inertial  effects  on  the  reaction 


dynamics  might  often  be  anticipated.  Especially  since  direct  experimental 
information  on  solvent  inertial  frequencies  are  lacking,  at  least  approxi¬ 
mately  reliable  analytic  formulations  of  u>o  are  clearly  desirable.  We  now 
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consider  specific  expressions  for  wD  ,  along  with  their  physical  basis. 

A  straightforward  formulation  of  v>0  has  been  given  by  Hynes,  who  wrote 

9 

(but  without  derivation)  the  following  expression 

“o  “  (2f»  +  ui  (2) 

gK 

In  Eq.  (2),  e,  and  ew  are  the  static  and  "infinite"  frequency  dielectric 
constants,  and  u{  is  the  "free  rotor"  frequency  which  for  linear  and 
spherical -top  molecules  is  given  by 

-  2^1/1  (2a) 

where  I  is  the  moment  of  inertia.  In  Appendix  A,  we  present  a  complete  de¬ 
rivation  of  Eq.  (2). 

The  presence  of  the  Kirkwood  "g"  factor  in  Eq.  (2),  gK ,  indicates  that 
solvent  structural  factors  are  important  in  determining  wo .  The  Kirkwood- 
Frolich  formula  is  an  adequate  means  of  calculating  gK  in  this  context^: 

<£.  -  0<2«,  +  O 


gK  “ 


+  2 V 


(3) 


where  N  is  the  Avogadro  Number  and  u  is  the  gas-phase  dipole  moment. 


Equation  (2)  reduces  to  an  expression  given  earlier  by  Calef  and  Wolynes 
for  polar,  non-polarizable  media  by  setting  gK  -  1  and  -  1  for  the  special 
case  of  a  rigid  rotor.  Their  equation  (72)  in  our  nomenclature  reads: 


Several  misprints  in  ref.  5  should  be  noted  in  this  connection.  Equation 
73  in  ref.  5  should  read  (in  their  notation): 


r  hr  <-jr>H 


Equation  77  should  read: 


r  - 


r  -1 

rR°T  u  ± 

-  (2«t  t)*5  exp(-/9AF* ) 

and  the  equation  for  k  in  Table  I  should  read: 

1 


k  *  —  'I +  h1  +  ^ 


K 


i  <.  C  T 


T-T-lV’ 


exp(-0AF*) 

2nC 


5 


V 


W°  “  (‘  I 


6y  ct 

7— r)li 


where  y,  the  reduced  concentration  of  dipoles  is  given  by 
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y  - 


Pf1  P 


(4) 


(4a) 


where  p  is  the  solvent  density.  Equation  (4)  can  be  rearranged  to  yield  a 
formula  having  the  same  structure  as  Eq.  (2): 


( 


21S  T 


-) 


(2c.  +  1) 


I  '  3gK 

Equation  (5)  can  also  be  expressed  in  the  following  form:' 

21S  T 


(5) 


wo  "  t<‘  x 


-)  7]»* 


(6) 


where 

eg  -  1  4w  pp20  1 

7  3y  31sT(e.  -  1)  "  1  +  g  ho  (6a) 

and  (!  is  a  normalizing  volume  element  and  hQ  is  the  pair  correlation  func- 
12 

tion.  This  formulation  reemphasizes  how  the  solvent  structure  (via  7) 
influences  wQ  and  hence  the  TST  reaction  rate.  The  use  of  Eq.  (5)  typically 
yields  ca  twofold  lower  estimates  of  «o  compared  with  Eq.  (2)  owing  to  the 
neglect  of  €„  in  the  former. 

It  is  interesting  to  note  that  the  expression  for  wo  given  in  Eq.  2  is 

also  similar  to  that  for  the  "dipolaron  frequency"  identified  by  Madden  and 
13 

Kivelson.  These  authors  developed  a  "three  variable  theory"  to  describe 
the  complex,  frequency- dependent  dielectric  constant,  e(u).  The  variables, 
Og  ,  flT  and  T  are  used  to  characterize  three  distinctly  different  processes 
that  occur  in  dipolar  fluids.  Two  variables,  fljj  and  flT ,  describe  dipolar 
rotational  and  librational  frequencies,  respectively;  the  last  variable,  T, 
describes  how  and  are  damped.  All  three  of  these  quantities  are  con- 
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sidered  in  the  context  of  a  time  correlation  function,  *(t) ,  which  is  related 

A 

to  «(w)  using  the  one-sided  Fourier- Laplace  transform,  JCiw: 

£(w>  -  €»  +  (<.  -  *.)  [1  +  Xiw  (*(t))]  (7) 

Although  all  three  variables  combine  to  dictate  the  polarization  dynamics 

A 

through  £ (w) ,  fij^  is  most  relevant  to  this  discussion.  This  quantity  repre¬ 
sents  an  inertial  frequency  (that  is,  a  frequency  which  is  inversely  propor¬ 
tional  to  the  moment  of  inertia)  that  characterizes  collective  dipolar  motions 
contained  in  the  short-time  behavior  of  5>(t).  For  linear  and  spherical  top 
molecules : 


<2«.  +  O 


2kjT 

I 


The  dipolaron  frequency  is  the  longitudinal  counterpart  of 


C i) 


2 

o 


(8) 


(9) 


as  can  be  seen  by  comparing  Eqs .  (8)  and  (9)  with  Eq.  (2)  (see  Appendix  A). 

Poley  originally  predicted  the  presence  of  a  spectral  feature  on  the 

basis  of  the  differences  between  and  the  square  of  the  refractive  index 

at  infrared  frequencies  for  "Debye"  flu: is  such  as  nitrobenzene.^  He  noted 

that  this  difference  was  proportional  to  n2  and  considered  this  unambiguous 

evidence  for  a  dielectric  loss  associated  with  molecular  reorientation. 

This  "Poley"  absorbance  has  since  been  observed  in  most  polar  liquids  that 

have  been  studied  in  the  far  infrared  region. Unfortunately,  u>o  cannot 

be  related  directly  to  the  frequency  of  the  Poley  band:  the  absorbance  is 

associated  with  the  librational  reorientation  of  solvent  dipoles  in  the  "cage" 

formed  by  neighboring  molecules,  which  is  a  transverse  (fiT ) ,  not  a  longi- 

13 

tudinal,  oscillation.  The  Poley  band  is  noted  here  since  it  constitutes 
one  of  several  high  frequency  dipole  reorientations  that  have  yet  to  be 
considered  in  the  context  of  their  relevance  (or  lack  thereof)  for  electron- 
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transfer  processes. 

The  practical  applicability  of  Eq.  (2)  is  subject  to  several  important 
constraints.  One  underlying  assumption,  of  course,  is  the  dielectric-con¬ 
tinuum  approximation.  Although  the  presence  of  molecular  moments  of  inertia 
and  dipole  moments  in  Eq.  (2)  implies  some  innate  consideration  of  solvent 
molecularity ,  the  many  countervailing  assumptions  that  are  made  to  assure 
tractability  may  lead  to  some  serious  uncertainties  in  the  numerical 
predictions.^  Nonetheless,  one  can  use  Eq.  (2)  as  well  as  molecular- dynamics 
simulations  to  provide  reasonable  estimates  of 

Estimates  of  obtained  from  Eq.  (2)  for  some  common  polar  solvents, 
encompassing  a  typical  range  of  dynamical  properties  (at  25°C) ,  are  listed 
in  Table  I  alongside  corresponding  values  of  r'1.  The  first  five  solvents 
each  exhibit  roughly  Debye  behavior  (i.e.  exhibit  only  a  single  clearcut 
dispersion  in  the  dielectric  loss  spectrum).  While  the  r"1  values  for  these 
solvents  vary  by  up  to  ca  AO  fold,  the  inertial  frequencies  are  seen  to  vary 
to  a  markedly  smaller  extent.  Given  that  with  the  exception  of  acetonitrile, 
tl1  <<  w0  •  one  expect  that  solvent  inertial  effects  upon  electron- 

transfer  dynamics  may  commonly  be  unimportant  at  ambient  temperatures.  How¬ 
ever,  acetonitrile  is  but  one  member  of  a  class  of  low-friction  solvents 
(also  acetone,  nitromethane ,  etc.),  for  which  inertial  effects  may  well  play 
a  significant  role. 

Corresponding  molecular  dynamics  simulations  of  are  as  yet  uncommon. 

Nevertheless,  such  calculations  yield,  for  example,  values  of  95  ps~ 1  for 

19  20 

water  and  23.5  ps'1  for  acetonitrile.  These  values  are  within  ca  twofold 
of  (albeit  sl’ghtly  higher  than)  the  corresponding  estimates  in  Table  I. 
Given  the  uncertainties  contained  in  both  Eq.  (2)  and  the  simulations,  this 
agreement  is  acceptable . 

For  larger  solvent  molecules,  however,  molecular  asymmetry  and  break¬ 
down  of  dielectric  continuum  theory  may  prove  a  more  serious  limitation  to 


8 


the  applicability  of  Eq.  (2).  Furthermore,  the  use  of  Eq.  (2)  together  with 

a  friction  kernel  for  the  rate  expression  that  does  not  include  the  short 

time  dynamics  of  the  longitudinal  polarisation  modes  brings  a  logical  incon- 

18 

sistency  into  the  calculations.  The  comple:  nature  of  twese  high-frequency 

16  21 

processes  has  been  described  in  detail.  ’  Further  insight  into  dielectric 
relaxation  processes  that  can  be  gleaned  from  inertial  effects  is  detailed 
in  Appendix  B  of  ref.  30.  For  the  calculations  that  follow,  we  choos  uo/2n 
values  o vei  the  range  5  to  25  ps"1,  reflecting  the  anticipated  range  for 
typical  solvents  (Table  I). 


2 .  Classical  Ad.abatic  Rate  expressions 

In  the  absence  of  solvent  friction  (i.e.  n  the  TST  limit)  the  classical 
adiabatic  barrier -crossing  frequency,  i/n ,  (i.e.  when  Tn  -  1,  x#1  -  1)  will 
equal  u>o/2»r  irrespective  of  the  shape  of  the  barrier  top  as  influenced  by 
the  electronic  coupling  matrix  element,  Hlz .  In  the  presence  of  solvent 
friction,  however,  »/n  should  generally  be  sensitive  to  the  barrier-top 

i  8 

curvature.-  Although  such  treatments  for  electron- transfer  reactions  have 
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been  discussed  previously,  ’  the  relevant  relationships  and  the  underlying 
phyp’cs  are  now  summarized  briefly  since  they  are  utilized  in  the  enduing 
numerical  calculations. 

In  general,  the  friction  at  the  barrier  top  is  probed  at  different  cros¬ 
sing  frequencies  than  the  function  experienced  by  a  particle  climbing  up 
(or  sliding  down)  the  well.  It  is  useful  to  make  a  piecemeal  parabolic  ap¬ 
proximation  to  the  true  potential -energy  surface:  fitting  two  parabolae  to 

q 

the  wells  and  an  inverted  parabola  to  the  barrier  top.  Since  u>  is  the 
characteristic  frequency  of  oscillations  in  the  wells,  simple  algebraic  ar¬ 


guments  yield  an  equivalent  frequency  near  the  barrier  top  given  by 
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2AG1 


-  1) 


*12 


(10) 
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where  AG*  is  the  "cusp"  barrier  height,  i.e.  in  the  limit  where  H12  -*  0. 

One  may  then  consider  the  overall  adiabatic  rate  constant  kA  as  a  composite 
of  the  rate  for  crossing  the  barrier  top,  1^,  and  that  for  diffusion  in  the 
wells,  kw: 


^  +  0-5  K1 


(ID 


Equation  (11)  is  a  restatement  of  the  Northrup-Hynes  "Stable  States  Picture" 

*  9  22 

of  barrier  crossing  dynamics,  applicable  if  AG  >  2-3  kjT.  ’ 

An  expression  for  in  the  classical  limit  is  the  Grote-Hynes  form- 

.  9,22,23 

ula: 


K  - 


2*  «b 


exp(  -AG*/^1) 


(12) 


where  \R  is  a  reactive  (or,  equivalently,  memory  renormalized)  frequency  that 
corresponds  to  the  largest  positive  root  of  the  equation 


w: 


"  AR  +  t}(\  ) 


(12a) 


given  the  frequency -dependent  friction,  »?(AR).  For  simplicity,  the  calcu¬ 
lations  here  consider  only  Debye  solvents,  so  that  the  friction  is  frequency 

,  18 

independent  and  »)  ■  or  r  L . 

The  treatment  of  kw  requires  a  different  approach.  Hynes  derived  an  in¬ 
tegral  expression  for  the  temporal  characteristics  of  a  Gaussian  (but  not 

9 

necessarily  Markov)  process  leading  to  kK.  In  the  Debye  limit  this 

5  9 

simplifies  approximately  to  ’ 

AG* 

k"  "  2rL1  (  «xp(-AG*/kBT)  (13) 


In  the  limiting  case  where  rL  -  0,  r\  -  0,  the  ratio  of  Ar/wb  appearing  in 
Eq.  (17)  is  unity.  As  then  kw  »  1^,  the  resulting  expression  for  kA  is  the 


usual  TST  result: 


10 


kA  *  k,,  -  («0/2ir)  exp  (  - AG*/kjT )  (14) 

More  generally,  in  the  presence  of  friction  kA  can  be  obtained  from  Eqs .  (11)- 
(13),  which  also  enables  the  classical  adiabatic  barrier-crossing  frequency 
[Eq.  (1)]  to  be  found  from 

‘'n  ”  kA  exp(AG*/kBT)  (15) 

III.  QUANTUM  CORRECTIONS  TO  THE  CLASSICAL  ADIABATIC  RATE  CONSTANT 

As  already  noted,  the  influence  of  the  nuclear  and  electronic  quantum 
corrections,  as  contained  in  Tn  and  k  x  ,  upon  the  reaction  rate  will  systemat¬ 
ically  be  greater  for  larger  u  values.  Unlike  i/  itself,  the  magnitude  of 
both  Tn  and  k  ^  in  the  absence  of  solvent  friction  is  expected  to  be  sensitive 
to  the  shape  of  the  barrier  top  as  well  as  the  wells.  Thus  diminishing  the 
extent  of  donor-acceptor  electronic  coupling  (i.e.  decreasing  H12)  will  en¬ 
large  at  a  given  ,  which  will  enhance  the  extent  of  nuclear  tunneling 
(i.e.  increase  Tn )  as  well  as  lessen  the  efficiency  of  electron  tunneling 
(i.e.  decrease  /c#l).  The  diminution  of  vn  for  a  given  barrier  shape  brought 
about  by  the  presence  of  solvent  friction  will  tend  to  decrease  the  magnitude 
of  these  quantum  corrections  (i.e.  depress  Tn ,  enlarge  *tl). 

Consequently,  then,  the  manner  in  which  such  quantum  corrections  can 
influence  solvent  dynamical  effects  in  electron  transfer  is  a  complex  inter¬ 
play  between  several  factors.  As  a  prologue  to  evaluating  typical  numerical 
consequences  of  such  quantum  effects,  we  now  outline  and  evaluate  some  can¬ 
didate  analytical  expressions. 

1.  Nuclear -Tunneling  Corrections 

A  convenient  expression  for  nuclear  tunneling  through  a  parabolic  barrier 
top  in  the  TST  limit  (i.e.  when  vn  —  u>o/2n)  is:^’^ 

wb  sinh  (A  uo/2\T) 
r"  “  V  sin  (A  wb/2kBT> 


(16) 
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At  higher  temperatures  (ftw  /2kjT  «1),  Eq.  (16)  can  be  approximated  by 


26 


(ft  wb/2kBT)  <*  Tk/2T) 

rn  *  sin  (ft  w^kjT)  =  sin  (*  Tk/2T)  (17) 

here  Tk  is  a  "characteristic  temperature"  of  the  system,  such  that  at  T  * 

Tk  the  probability  for  transitions  above  the  barrier  and  those  through  it 

are  equal.  This  latter  formulation  is  of  particular  value  as  a  diagnostic 

tool:  if  T  <  |  Tk,  large  nuclear  tunneling  effects  are  predicted  whereas 

26 

for  T  >  2Tk  tunneling  effects  will  be  small  or  negligible. 

Both  formulas  given  above  face  restrictions  concerning  the  admissible 

degree  of  barrier  curvature;  specifically,  they  require  that  hu>b  <  2>rkBT 

to  assure  that  Tn  remains  a  positive  definite  quantity.  These  approximate 

27 

formulations  also  face  difficulties  when  Tn  is  large  (say,  >  5-10).  Such 
a  circumstance  indicates  that  tunneling  from  deep  within  the  well  is  sig¬ 
nificant,  complicating  the  analysis. 

An  expression  for  the  tunneling  correction  that  is  appropriate  for  cusp- 

28 

like  barriers  (i.e.  when  H12  is  small)  has  been  given  by  Holstein: 


r 


n 


sinh  (ft  a>o/2kgT)  4AG*  ft  vo  ft  u>o 

(ft  we/2kBT)  eXpi'(ft“^)  [tanh(_4k^f)' 


(18) 


Equation  (18)  should  be  valid  over  a  wide  temperature  range. 

Introducing  friction  into  the  nuclear- tunneling  correction,  as  is  needed 

here,  requires  a  rather  more  involved  analysis.  A  lucid  discussion  of  the 

quantum  corrections  to  the  transition  state  theory  and  the  generalized  quantum 

24 

Kramer's  theory  has  been  provided  by  Hanggi.  As  noted  there,  and  by 
29 

Wolynes,  the  ratio  of  AR/u>b  that  distinguishes  the  classical  adiabatic 


rate  constant  from  [Eq.  (12)]  also  serves  to  renormalize  the  char¬ 
acteristic  temperature,  Tk,  used  in  the  calculation  of  Tn  [Eq.  (17)].  Stated 
differently,  by  regarding  AR  as  a  memory  renormalized  frequency,  the  presence 


12 


of  energy  dissipation  acts  to  change  the  path  and  frequency  of  barrier  cros- 
24 

sing.  In  the  present  analysis,  the  nuclear- tunneling  correction,  r  ,  is 
taken  into  account  entirely  in  the  "barrier- top"  component,  ,  of  the  overall 
rate  constant  kA  [Eqs.  (11),  (12)]: 

K  -  fn  (AR/Ub)(Wo/2*)  exp(-4G*ABT)  (19) 

The  overall  nuclear- tunneling  correction,  Tn,  in  Eq.  (1)  is  thereby  smaller 
than  f  ,  being  weighted  by  the  relative  magnitudes  of  k^  and  kw: 


<2^  +  k„)  Fn 
2rn  ^  +  kw 


(20) 


While  this  formulation  satisfactorily  considers  the  nuclear- tunneling  cor¬ 
rection  to  kjj,  it  can  also  be  questionable  if  Fn  is  large  since  a  tunneling 
correction  to  kw  would  also  then  be  necessary. 

29 

A  useful,  albeit  cumbersome,  formulation  of  Tn  was  derived  by  Wolynes 

25  24 

and  by  Daknovskii  and  Ovchinnikov  (and  recapitulated  lucidly  by  Hanggi  ) : 


(21) 


(21a) 


(21b) 


given  6  ■  and  r)  ■  The  variables  A+  and  A*  are  familiar  as 

the  two  poles  obtained  from  the  Grote-Hynes  analysis  for  the  case  of  freq- 

9  22  23 

uency- independent  friction.  ’  ’  Equation  (21b)  establishes  the  (memory 

renormalized)  frequencies  characterizing  fluctuations  about  the  minimum. 


Note  that  Eq.  (21)  has  simple  poles  at  A+/0  -  1  and  at  A ~/6  -  1.  Moreover, 

Eq.  (21b)  becomes  complex  for  ij  <  2wo . 

It  can  be  shown  that  Eq.  (21)  reduces  analytically  to  Eq.  (16)  in  the 
30 

zero-friction  limit.  Given  that  Eq.  (16)  can  be  approximated  under  s^me 
conditions  by  Eq.  (17),  this  suggests  that  Eq.  (21)  might  also  be  approximated 
by  the  memory  renormalized  version  of  Eq.  (17): 

Ar  sinh(ft  w0/2kBT) 

? n  “  ^7  sin  (ft  V2V>  (22) 

The  effect  of  friction  upon  Tn  is  therefore  considered  to  arise  simply  from 

the  diminution  of  the  net  velocity  along  the  reaction  coordinate  close  to 

the  barrier  top  from  the  TST  value,  u>b  ,  to  AR  .  This  treatment  of  friction 

effects  in  terms  of  renormalized  velocities  is  closely  analogous  to  procedures 
commonly  used  to  describe  such  effects  upon  reaction  nonadiabaticity  (vide 
infra).7  We  have  found  Eq.  (22)  to  be  a  good  approximation  to  Eq.  (21)  (with¬ 
in  ca.  20%)  over  the  range  of  parameters  encountered  in  the  numerical  calcu¬ 
lations  described  below.  Equation  (22)  has  the  virtue  of  allowing  the  exam¬ 
ination  of  the  complete  rate-friction  regime,  using  a  single,  physically 

transparent,  formula.  Another  simple  expression  which  provides  a  reasonable 

9 

approximation  to  Eq.  (21)  is  (cf  ref.  31): 

fn  »  exp((ft/kBT)2(ujg  +  u2)/24)  (23) 

In  view  of  Eq.  (10),  Eq.  (23)  can  be  rewritten  as 

Fn  =  exp((ft/kBT)2AG*  w2/12H12)  (24) 

Equation  (24)  emphasizes  the  sensitivity  of  fn  to  wo,  H12  and  the  temper¬ 
ature  . 

Appendix  B  presents  a  brief  numerical  comparison  of  nuclear-tunneling 
factors  extracted  using  some  formulae  presented  in  this  section. 
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2 .  Nonadiabaticity  Corrections 

The  anticipated  influence  of  reaction  nonadiabaticity  [i.e.  <  1  in 

Eq.  (1)]  on  the  ET  rate-solvent  friction  dependence  is  well  docu- 
1  7  9  32  33 

mented.  '  ’  ’  ’  The  prediction  that  the  extent  of  this  dependence  can 

be  sensitive  to  both  the  donor-acceptor  electronic  coupling  and  the  magnitude 
of  the  friction  has  recently  received  direct  experimental  verification. ^ 
Since  we  have  discussed  the  treatment  of  such  electronic  nonadiabaticity 
effects  in  detail  elsewhere,^  our  remarks  here  are  confined  to  a  brief  summary 
of  simplified  formulae  relevant  to  the  numerical  calculations  described  below . 

A  generally  useful  expression  for  nonadiabatic  effects  upon  the  unimolec- 
ular  ET  rate  constant  in  the  presence  of  solvent  friction  is^ 

«.l  k?ST 

U.t -  (25) 

1  -  Si  +  KikTSrAA] 

where  kTST  is  the  rate  constant  corresponding  to  the  adiabatic  TST  limit, 
and  K'j'  is  related  to  the  overall  electronic  transmission  coefficient  1  by 

*,l  “  +  «.i<l  -  *„)]  (26) 

where  xj|  -  kA/kTET  .  The  k#1  values  are  related  to  the  extent  of  electronic 
coupling  by^ 


*#1  “  4»r3/27T/(+  A7r3/27l)  (27) 

where 

7T  -  iH^IVAftu^AG-kjT)1*  (27a) 

A  particularly  straightforward,  albeit  more  approximate,  formulation 

of  *gl  which  is  applicable  either  in  the  presence  or  absence  of  friction  is 

7  32 

the  Landau-Zener  multiple-crossings  expression  ’ 

2  [  1  -  exp( -v %  x  /2*/„  ) ) 

*#1  2  -  exp( j /2un ) 

where  the  "electronic  frequency  factor"  i/#1  is  given  by 


(28) 


(28a) 


-  K\z(*/U>hG* h2\T)h 
Since  the  nuclear  frequency  factor  describes  the  net  velocity  for  classical 
adiabatic  passage  over  the  barrier  top,  it  accounts  simply  for  the  effect 
of  solvent  friction  upon  <c#1  [cf  Eq.  (22)  for  nuclear  tunneling].  Equation 
(28)  therefore  provides  a  useful  means  of  interpolating  between  the  adiabatic 
limit  (where  -+  1)  to  the  nonadiabatic  limit,  where  a  H22. 

IV.  NUMERICAL  RESULTS 

Having  outlined  reasonably  versatile  means  of  formulating  the  influence 
of  nuclear  and  electronic  quantum  corrections  on  the  ET  rate  both  in  the 
presence  and  absence  of  solvent  friction,  it  is  of  interest  to  examine  the 
combined  numerical  consequences  of  such  effects,  especially  regarding  their 
influence  upon  the  anticipated  rate- friction  dependence  in  the  vicinity  of 
the  inertial  limit. 

Figures  1  and  2  show  how  the  friction  influences  the  magnitude  of  the 
nuclear  tunneling  correction,  F  ,  as  obtained  from  Eq.  (21),  with 
uQ  -  5  x  1012  or  2.5  x  1013  s'1,  respectively,  T  -  298  K,  and  AG*  -  5.0  kcal 
mol'1.  The  horizontal  dashed  lines  represent  the  TST  limiting  values  of  fn 
for  the  electronic  matrix  coupling  (H12)  values  indicated  [as  given  by  Eq. 
(16)],  while  the  corresponding  solid  curves  include  the  effects  of  solvent 
friction.  Comparing  the  ?n  values  given  in  Figs.  1  and  2  illustrates 
the  sensitivity  of  the  nuclear  tunneling  correction  to  the  magnitude  of  «o 
as  well  as  to  ^  (i.e.  to  H12),  particularly  in  the  low  friction  (high  r^1 ) 
regime.  This  sensitivity  was  Initially  surprising  to  us.  However,  the 
observation  becomes  readily  understandable  upon  examining  the  simplest  approx¬ 
imate  formula  for  Fn,  Eq.  (24),  in  that  this  relation  exhibits  an  exponential 
dependence  of  Tn  upon  H12 .  The  inclusion  of  friction  effects,  however,  re¬ 
duces  the  magnitude  of  fn  as  well  as  its  dependence  upon  H12. 

In  a  similar  fashion,  Fig.  3  illustrates  how  the  presence  of  solvent 
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friction  can  diminish  the  dependence  of  Fn  upon  the  reciprocal  temperature 
(1C1)  as  calculated  using  Eqs.  (16)  and  (21).  The  dotted  curve  representing 
the  nuclear -tunneling  correction  in  the  TST  limit  (for  H12  -  0.2  kcal 
mol  *1)  is  seen  to  be  significantly  more  temperature  depend  nt  than  are  the 
corresponding  solid  curves  which  illustrate  the  influence  of  friction.  The 
latter  show  that  Tn  decreases  markedly  with  increasing  friction;  the  three 
solid  curves  refer  to  rLwo  values  of  2,  5,  and  20,  being  fixed  at  2  x  1013 
s'1.  The  dashed  curve  in  the  upper  left-hand  corner  of  Fig.  3  illustrates 
further  the  greater  importance  of  nuclear  tunneling  for  more  cusp-like  bar¬ 
riers;  here  H12  -  0.05  kcal  mol'1  with  rLa>0  -  20.  The  TST  results  shown  in 
Fig.  3  merit  a  note  of  caution,  however:  Eq.  (16)  becomes  a  progressively 
worse  approximation  as  the  temperature  approaches  the  divergence  point  at 
T  —  fto;b/2jrkB.  This  limitation  arises  from  the  presumption  that  the  barrier 
top  can  be  approximated  by  an  inverted  parabola. 

Figures  4-7  illustrate  how  the  combination  of  nuclear  and  electronic 
quantum  corrections  can  affect  the  dependence  of  the  reaction  rate  on  solvent 
friction,  in  the  form  of  logarithmic  plots  of  k#t  versus  t^1.  (Although  we 
have  previously  evaluated  bimolecular  rate  constants  by  spatial  integration 
of  ktt  values  so  to  explore  detailed  friction- dependent  effects^,  attention 
is  restricted  here  to  the  latter  unimolecular  quantity  for  simplic  ity)  . 
The  results  in  Fig.  4  were  calculated  by  taking  H12  -  0.2  kcal  mol'1,  AG*  - 
5.0  kcal  mol'1,  and  T  -  298  K  (see  the  Figure  Caption  for  calculational  de¬ 
tails).  The  inertial  frequency  was  also  fixed  at  5  x  1012  s'1;  this  value 

13 

was  selected  as  being  a  typical  value  for  polar  solvents  (vide  supra). 
The  range  of  rL  values,  1  x  1011  s' 1  to  5  x  1012  s' 1 ,  reflects  that  commonly 
encountered  in  polar  Debye-like  media. ®  The  solid  log  k#t  -  log  r^1  curve 
was  obtained  without  quantum  corrections;  the  curvature  seen  towards  larger 
r'Lx  arises  entirely  from  the  progressively  larger  influence  of  the  TST  limit 
set  by  the  chosen  inertial  frequency.  The  dotted  curve  includes  the  influence 


17 


of  nuclear  tunneling,  and  the  dashed  curve  the  combined  effects  of  nuclear 
tunneling  and  electronic  nouadiabaticity .  (Admittedly,  in  real  experimental 
systems  alterations  in  rL  brought  about  by  solvent  substitution  also  yield 
significant  variations  in  AG*  and  wQ  .  Nevertheless,  holding  the  latter  para¬ 
meters  fixed,  as  in  Figs.  4-7,  enables  their  influence  on  the  rate-friction 
functionality  to  be  more  readily  discerned) . 

The  accelerating  influence  of  the  former  quantum  effect  is  seen  to  be 
small  under  these  conditions,  being  swamped  by  the  progressive  rate  retarda¬ 
tion  seen  towards  larger  r'1  that  arises  from  reaction  nonadiabaticity .  The 
latter  effect  has  been  discussed  in  detail  recently.^ Not  surprisingly, 
the  extent  of  friction- dependent  rate  retardation  due  to  reaction  nonadia¬ 
baticity  increases  markedly  as  Hl2  decreases.  While  the  importance  of  nuclear 
tunneling  also  increases  under  these  conditions,  reaction  nonadiabaticity 
tends  to  constitute  the  dominant  quantum  effect  for  a  wide  range  of  H12  values 
for  the  barrier  parameters  as  in  Fig.  4. 

Significantly  different  behavior,  however,  can  be  obtained  under  other¬ 
wise  similar  conditions  but  by  choosing  somewhat  larger  u>0  values.  This  is 
illustrated  in  Fig.  5,  which  contains  log  k#t  -  log  r'1  plots  obtained  as 
in  Fig.  4  but  for  a  fivefold  larger  w0  value,  2.5  x  1013  s'1.  This  inertial 
frequency  is  characteristic  of  relatively  small  solvent  molecules  such  as 
acetonitrile  (vide  supra).  Inspection  of  Fig.  5  shows  that  while  the  deceler¬ 
ating  effect  of  reaction  adiabaticity  still  tends  to  be  greater  than  the 
accelerating  effect  of  nuclear  tunneling  (compare  solid,  dashed  curves), 
the  magnitude  of  the  latter  is  now  comparable  to  the  former  (compare  solid, 
dotted  curves).  This  is  especially  the  case  for  high  or  moderate  solvent 
friction  (i.e.  lower  r'1  values). 

Given  that  the  numerical  importance  of  nuclear  tunneling  is  sensitive 
to  the  temperature  (Fig.  3),  it  is  of  interest  to  examine  such  combined  quan¬ 
tum  effects  on  the  rate- friction  behavior  at  lower  temperatures.  For  this 


purpose,  Figs.  6  and  7  contain  log  k#t  -  log  t'1  plots  calculated  as  in  Figs. 
4  and  5,  respectively,  but  with  T  -  150  K.  The  markedly  enhanced  influence 
of  nuclear  tunneling  is  clearly  evident  in  these  results.  In  particular, 
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Fig.  7  (for  which  u>o  -  2.5  x  1013  s'1)  shows  that  the  rate -decelerating  non- 
adiabatic  effects  are  largely  overshadowed  by  nucleai  tunneling,  yielding 
k#t  values  (dashed  curve)  that  are  larger  than  the  corresponding  classical 
quantities  (solid  curve)  throughout  the  range  of  solvent  friction  shown. 

General  Conclusions 

One  overall  conclusion  which  emerges  from  such  numerical  calculations 
is  that  the  inclusion  of  such  quantum  effects  can  mask  partially  or  even 
completely  the  influence  of  the  solvent  inertial  TST  limit  on  the  reaction 
rate.  Electron  tunneling  rather  than  nuclear  dynamics  is  often  expected  to 
at  least  partly  control  the  net  barrier-crossing  velocity  [i.e.  /ctl»/n  -  i/el 
rather  than  «el^n  -  va  ,  Eq.  (28)]  once  the  latter  approaches  the  anticipated 
inertial  limit  frequencies,  ca  1013  s'1.  Nevertheless,  some  influence  of 
the  inertial  dynamics  can  be  anticipated  for  systems  displaying  moderate  or 
strong  donor-acceptor  electronic  coupling,  say  H12  >  0.2  kcal  mol'1.  While 
the  compensatory  influence  of  nuclear  tunneling  associated  with  solvent  iner¬ 
tia  is  typically  milder,  this  factor  can  mask  the  rate-limiting  effects  of 
both  reaction  nonadiabaticity  and  the  TST  inertial  frequency.  Given  that  the 
quantum  effects  associated  with  nuclear  tunneling  and  reaction  nonadiabaticity 
arc  inherently  compensatory  and  friction-dependent,  the  common  experimental 
stracegy  of  examining  the  ET  rate-solvent  friction  dependence^  may  not  provide 
a  clear  indicator  of  TST  limiting  effects.  Indeed,  the  occurrence  of  the 
classical  solvent -controlled  TST  limit,  whereby  the  barrier-crossing  freq¬ 
uency  equals  u>Q/2n,  should  require  the  presence  of  substantial  elec¬ 
tronic  coupling  (i.e.  large  H12)  so  to  both  maintain  reaction  adiabaticity 


and  diminish  the  importance  of  nuclear  tunneling. 

While  nuclear  tunneling  involving  solvent  inertial  motions  may  often 
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be  small  or  negligible  at  room  temperature,  Fig.  7  illustrates  the  potential 
importance  of  this  phenomenon  at  lower  temperatures.  Electron- transfer 
kinetic  experiments  of  this  type  concerned  with  solvent  dynamical  phenomena 
have  apparently  yet  to  be  performed.  Nevertheless,  there  are  several  prac¬ 
tical  reasons  to  earmark  such  studies  for  future  consideration,  not  the  least 
of  which  is  the  markedly  diminished  reaction  rates  expected  under  these  condi¬ 
tions  . 

Quite  apart  from  nuclear  tunneling  involving  solvent  rotational/libra- 
tional  modes,  substantial  quantal  effects  involving  high-frequency  reactant 
vibrational  and  related  solvent  modes  can  often  be  anticipated.  For  suf¬ 
ficiently  high  frequencies  (say  >  5  x  1013  s'1),  especially  for  cusp-like 
barriers  (i.e.  small  H12),  nuclear  tunneling  involving  such  modes  can  be 
sufficiently  important  even  at  ambient  temperatures  so  that  their  effective 
contribution  to  the  activation  barrier  is  muted  or  even  essentially  removed. 
Interesting  recent  calculational  examples  of  such  phenomena  include 
Fe(OH2)3+/2+  self  exchange  in  aqueous  media,  for  which  the  presence  of  sub¬ 
stantial  nuclear  tunneling  was  deduced  for  hydrogen  atoms  motions  involving 

34 

the  solvent  as  well  as  the  aquo  ligands.  Along  with  the  high  wo  values 
estimated  for  water  (-101*  s'1,  vide  supra),  sizable  nuclear- tunneling  effects 
might  often  be  expected  in  this  solvent.  Such  effects  will  be  lessened  sub¬ 
stantially,  however,  for  processes  featuring  strong  donor -acceptor  orbital 
overlap,  such  as  the  cobaltocenium-cobaltocene  self -exchange  reaction  recently 
evaluated  in  water  as  well  as  other  solvents  in  our  laboratory .  ®  ^  ^  More¬ 
over,  the  reaction  dynamics  may  well  be  largely  overdamped  in  water  (i.e.  sol¬ 
vent  friction  is  prelevant),  since  r'1  «  u>0  (Table  I). 

The  uncertainties  in  estimating  the  solvent  inertial  frequency,  noted 
above,  could  be  construed  as  providing  a  major  impediment  to  numerical  cal¬ 
culations,  especially  involving  larger  and  more  asymmetric  molecules.  Fortun¬ 
ately,  however,  the  commonly  anticipated  prevalence  of  overdamped  motion 
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18 

for  such  solvents  makes  such  uncertainties  less  serious  in  practice.  Of 
greater  anticipated  significance  is  the  occurrence  of  additional  dissipa¬ 
tive  relaxations  at  higher  frequencies  than  r^1.  These  oft-present  com¬ 
ponents  can  diminish  greatly  the  effective  solvent  friction,  yielding  substan¬ 
tial  rate  accelerations  beyond  that  anticipated  for  frequency- independent 
9  18  35  36 

friction.  ’  ’  ’  A  quantitative  understanding  of  the  influence  (or  lack 

thereof)  of  such  high-frequency  frictional  components,  however,  will  inevit¬ 
ably  require  detailed  consideration  of  quantum  effects. 
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Appendix  A;  Derivation  of  Ea .  (2)  and  Related  Considerations 

Equation  (2)  can  be  generalized  to  include  a  somewhat  more  involved 

definition  of  «f  than  Eq.  (2a).  The  free  rotor  frequency,  wf ,  is  defined 
9,37 

as : 

w2  —  <u‘w>  -  <(u'u>)2>  -  <u2>  A.l 

where  u  is  the  molecular  angular  velocity  vector,  and  u  is  the  unit  vector 

that  lies  along  molecular  dipole  moment  vector,  y. .  Here,  we  will  explore 

9  37 

the  derivation  more  fully,  following  Hynes  and  Brot.  Evaluating 

u2  ■  <ww>  -  <(u-w)2>  term  by  term,  we  find 

W£  “  <wx  +  Wy  +  Ul  -  <UXWX  +  Uy‘Wy  +  U2'W2)Z>  A.  2 

Evaluating  the  square,  and  taking  advantage  of  the  fact  that  cross  terms  in 

37 

Eq .  (A. 2)  vanish,  allows  the  simplification: 

Ul  “  ^x  +  Wy  +  Wz  '  (<U2  +  Uy)Wx  +  'Uy  +  Ux)wx  +  <Uz  +  u^)Wy  ]>  A.  3 

Since  u  is  a  unit  vector,  u2  +  u2  +  u2  -  1,  where  ux,  uy  and  ur  are  the  pro¬ 
jections  of  the  unit  vector  along  the  dipole  moment  onto  the  principal  axes 

of  inertia,  the  formula  for  w2  in  equation  A.  3  (after  inserting  the  equ-parti- 

3  7 

tion  values  for  w,  i.e.  <*>2  -  kfiT/IQ)  becomes: 


(u2  + 


uy> 


Uz) 


(U2  +  U2) 


+ 


A.  4 


or,  restated  in  more  compact  form, 


,  M  *  »a 

•4  -  V  2  (—3 -  )  (i;1)  A. 5 

a  p 

where  nz  -  (T  connotes  the  vector  transpose)and  q  -  x,y,z  This  more 

general  expression  reduces  to  Eq.  (2a)  in  the  text  for  linear  and  spherical - 
top  molecules . 
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The  three  individual  components  of  the  complete  dipole  moment  vector 

can  be  uerived  from  the  Stark  effect  splitting  of  the  'gas  phase]  rotational 

spectra  for  the  molecule.  (Dipole  moments  are  often  markedly  different  in 
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condensed  phase  than  in  vacuum  so  the  use  of  these  parameters  generates  some 
uncertainty  in  estimating  uo . )  The  principal  moments  of  inertia  may  be  ex¬ 
tracted  from  tabulations  of  the  three  rotational  constants,  RQ  (usually  called 

37 

A,  B  and  C,  respectively),  using  the  relation: 


R 


a 


41Qc 


A. 6 


where  c  is  the  speed  of  light. 

The  dielectric  parameter  terms  in  Eq.  (2)  arise  from  the  form  of  the 
dielectric  response  function,  expressed  as  the  dielectric  permittivity 


E(s) 


37,39 


E(s)  -  {1+ 


3c. 


2c.  +  e„ 


Sk 

[^(-g(t)) 


l]}*1 


2cs  +  c(s) 

3e.  Sk 


i  (-g(t)) 


A.  7 


where  g(t)  -  dg(t)/3t.  Equation  A. 7  utilizes  the  notion  of  the  time -dependent 
Kirkwood  "g"  factor  -  a  collective  function  describing  the  structural  correla¬ 
tions  in  the  ensemble  of  solvent  molecules.  If  the  structural-correlation 
function  decays  exponentially  in  time,  /(-g(t))  a  1/(1  +  sr*) ,  in  which  case 

a 

E  becomes  a  simple  exponential  with  the  relaxation  time,  r*.  characterizing 

the  decay  of  the  structural  correlation  function  g(t);  r*  is  related  to  che 

37  39 

Debye  relaxation  time  rD  by:  ’ 


3c. 


TD  " 


2ct  +  cc 


A .  E 


As  pointed  out  by  Hynes,  g(t)  obeys  a  Langevin  equation  having  the  usual 
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form: 


g(t)  -  -J  dr  K(t  -  r)  g(r ) 


where  the  kernel,  K(t) ,  has  the  Laplace  transform: 


K(s)  - 


s  +  £l(s)  gK[(s  +  £l(s)] 


The  frictional  term  £,(s)  contains  the  dielectric  loss  information.  The 


.2  /_  9 


kernel,  K(0)  in  Eq.  (A.  10),  equals  w2/ gk .  Inserting  this  into  the  overall 

A 

expression  for  the  dielectric  response  function  E(s)  as  given  in  equation 
A.  7,  the  high  frequency  term  (multiplied  by  s2)  can  be  identified  as  an  "iner¬ 
tial  limited"  transverse  frequency  if  we  assume  a  friction  coefficient  having 
the  form  fT(s)  -  w2  rD(s).  Assuming  rD(s)  -  rD  (i.e.  Debye  behavior), 

A 

E(s)  can  then  be  reexpressed  as: 


E(s)  -  {1  +  +  sr  r 

wT 


where  wT  is  now  identified  as: 


2«,  + 

wt  “  3«g 


A. 11a 


One  then  extracts  the  corresponding  longitudinal  quantity,  u>  ,  using 


[cf  Eq.  (9)] : 


which  yields 


wo  "  wt  <e./e«) 


(2«,  +  tm)  u>\ 

““  Hnr - 

Eq.  (A. 13)  is  identical  with  Eq.  (2)  in  the  text.  It  is  also  worth  noting 
that  the  proportionality  coefficient  relating  the  transverse  and  longitudinal 
inertial  frequencies  wT  and  w0  in  Eq .  (A. 12),  ( e , / f ^ ,  is  similar  to  but 
not  identical  with  the  well-known  factor  («,/<„)  relating  the  transverse 
and  longitudinal  overdamped  relaxation  frequencies  r'1  and  t*1.  A  detailed 


derivation  of  Eq.  (A. 12)  along  with  a  discussion  of  the  underlying  physics 
responsible  for  these  differences  is  to  be  found  in  Appendix  A  of  ref.  30. 
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Appendix  B:  A  Numerical  Comparison  of  Nuclear -Tunneling  Corrections 

We  outlined  above  a  number  of  different  approximate  formulations  of 
the  nuclear  tunneling  correction  factor  Tn.  Table  B.I  presents  a  numerical 
comparison  of  Tn  values  using  Eqs.  (16),  (17),  and  (24)  for  three  different 
values  of  H1£  and  «o .  It  is  seen  that  all  three  expressions  yield  similar 
results  when  wQ  —  5  x  1012  s’1.  However,  significant  discrepancies  (even 
divergences!)  are  obtained  at  larger  values  of  wQ,  particularly  for  small 
values  of  H12. 


Table  B.I 


A  Comparison  of  Nuclear  Tunneling  Correction  Factors 
in  the  Transition  State  Theory  Limit  at  298  K 


Approximation 

Hia  3  x 
(kcal  mol  A) 

w  -  5xl012 
°  (s’1) 

w  -  lxlO13 
(s’1) 

u>  -  2.5x10 
(s’1) 

Equation  16 

0.05 

1.15 

1.87 

c 

0.125 

1.06 

1.34 

10.04 

0.2 

1.04 

1.26 

2.92 

Equation  17 

0.05 

1.15 

1.86 

c 

0.125 

1.06 

1.33 

9.87 

0.2 

1.03 

1.25 

2.87 

Equation  24 

0.05 

1.15 

1.79 

5.41 

0.125 

1.06 

1.32 

3.93 

0.2 

1.04 

1.25 

2.35 

a  The  cusp  barrier  height  was 

“ft 

taken  as  AG  - 

C 

5.0  kcal  mol' 

l 

®  The  barrier 

top  frequency,  wt ,  was  calculated  using  Eq. 

(10). 

c  No  numerical 

estimate  could 

be  obtained  owl 

ng  either  to 

divergence  in 

expression  or  violation  of  boundary  conditions  for  the  approximation. 
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TABLE  I  Comparison  Between  Inertial  Frequencies  Estimated  from  Eq.  (2)  and 
Inverse  Longitudinal  Relaxation  Times  for  Some  Common  Solvents  at 
25°C. 


Solvent 

ps'1 

b 

Wo 

ps'1 

Acetonitrile 

4 

11 

d2o 

1.9 

40 

Dimethylsulfoxide 

0.5 

9.5 

Benzonitrile 

0.2 

4 

Hexamethylphosphoramide 

0.11 

-4 

Methanol 

(0.135) 

11 

Ethanol 

(0.033) 

9.5 

a  Inverse  longitudinal  relaxation  times 

from  dielectric  loss  data  (see  ref. 

for  solvent  indicated,  as  obtained 

6  for  data  sources) .  Values  for 

methanol  and  ethanol  (given  in  parentheses)  refer  to  large - amplitude , 
longer  time,  portion  of  multicomponent  dielectric  disposition. 

Solvent  inertial  frequency,  estimated  from  Eq.  (2)  (see  ref.  18  for 
data  sources). 


b 
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FIGURE  CAPTIONS 
Fie.  1 

Nuclear- tunneling  correction  to  rate  constant  as  a  function  of  log  r ‘ 1  , 
where  rL  is  the  solvent  longitudinal  relaxation  time,  for  three  different 
electronic  matrix  coupling  elements  H12  as  indicated,  calculated  from  Eqs . 
(10)  and  (21) .  The  dashed  horizontal  lines  represent  the  TST  limiting  value 
(i.e.  for  rL  -*  0) ,  calculated  from  Eqs.  (10)  and  (16).  Other  parameters: 

-  5  x  1012  s’1,  AG*  -  5.0  kcal  mol'1,  T  -  298  K. 

EliU-2. 

As  in  Fig.  1,  but  for  H12  values  indicated  with  a  higher  inertial  limit¬ 
ing  frequency,  u>o  -  2.5  x  1013  s'1. 

m.,  3 

Nuclear- tunneling  correction  factors  as  a  function  of  the  reciprocal 
absolute  temperature.  Solid  curves  refer  to  rL  values  (in  order  of  decreasing 
Tn )  equal  to  2/wq  ,  5/u>0  ,  and  20/wo ,  with  H12  -  0.2  kcal  mol'1,  calculated 
by  using  Eqs.  (10)  and  (21).  Dashed  curve  refers  to  r,  -  20/w  and  H.  ,  - 
0.05  kcal  mol'1.  Dotted  curve  represents  the  TST  result  (i.e.  for  rL  -►  0) 
for  H12  -0.2  kcal  mol'1  obtained  from  Eqs.  (10)  and  (16).  Other  parameters 
wQ  -  2  x  1013  s'1,  AG*  -  5.0  kcal  mol'1. 

EifU-A 

Logarithmic  plot  of  the  unimolecular  ET  rate  constant,  k  ,  versus  the 
inverse  longitudinal  solvent  relaxation  time.  Solid  curve  is  the  classical 
adiabatic  result,  obtained  from  Eqs.  (1),  (10) -(13),  and  (15),  with 
<*>o  -  5  x  1012  s'1,  H12  -  0.2  kcal  mol'1,  AG*  -  5.0  kcal  mol'1,  T  -  298  K. 
Dotted  curve  includes  nuclear  tunneling  correction  as  calculated  from  Eq. 
(21).  Dashed  curve  includes  additionally  the  electronic  transmission  co¬ 
efficient  (i.e.  reaction  nonadiabaticity)  calculated  from  Eq.  (28). 
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As  in  Fig.  4,  but  for  «o  -  2 . 5  x  1013  s'1  and  T  -  150  K. 
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